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Pauli-spin blockade (PSB) is a transport phenomenon in double quantum dots that allows for a 
type of spin to charge conversion often used to probe fundamental physics such as spin relaxation 
and singlet-triplet coupling. In this paper we theoretically explore Pauli-spin blockade as a function 
of magnetic field B applied parallel to the substrate. In the well-studied low magnetic field regime, 
where PSB occurs in the forward (1,1) -A (0, 2) tunneling direction, we highlight some aspects of 
PSB that are not discussed in detail in existing literature, including the change in size of both bias 
triangles measured in the forward and reverse biasing directions as a function of B. At higher fields 
we predict a crossover to “reverse PSB” in which current is blockaded in the reverse direction due 
to the occupation of a spin singlet as opposed to the traditional triplet blockade that occurs at low 
fields. The onset of reverse PSB coincides with the development of a tail like feature in the measured 
bias triangles and occurs when the Zeeman energy of the polarized triplet equals the exchange energy 
in the (0,2) charge configuration. In Si quantum dots these fields are experimentally accessible; thus, 
this work suggests a way to probe singlet to triplet relaxation mechanisms in quantum dots when 
both electrons occupy the same quantum dot. 


Since its discovery, Pauli-spin blockade (PSB)[1] has 
been a valuable tool for probing fundamental physics. Its 
use in spin to charge conversion has led to investigations 
of spin Ti relaxation times [2, 3], electron spin couplings 
to lattice nuclear spins[4], and spin-orbit effects[5]. PSB 
has also received a lot of attention from the quantum 
information community due to its use in read-out and 
initialization of various electron spin states [6-8]. The 
majority of studies involving PSB have focused on the low 
magnetic field regime. In this manuscript we investigate 
the qualitative behavior of PSB at all magnetic fields. 
In doing so we identify experimentally accessible regimes 
where PSB has yet to be studied and new physics is likely 
to be found. 

In a typical double quantum dot (DQD) device two 
quantum dots are coupled to each other as well as 
to two electron reservoirs (see figure 1(a)). The dots 
are also capacitively coupled to one or more gate elec¬ 
trodes which can be used to control the dot poten¬ 
tials. The charge configuration of the system is char¬ 
acterized by a pair of numbers (n, m) corresponding to 
the number of electrons on the left and right dot, re¬ 
spectively. Transport through quantum dots can be un¬ 
derstood by considering the chemical potentials of the 
quantum dots, /i/(n, m) = E(n,m ) — E(n — 1, m) and 
fi r {n , m) = E(n, m) — E(n , m — 1), where E denotes the 
total energy of the system. When the appropriate chemi¬ 
cal potentials form a energetically downhill path between 
the Fermi levels of the two leads, Ull(rl)i electron tun¬ 
nelling through the device is energetically favorable and 
leads to current flow[9]. Near the (1,1) «a (0,2) charge 
transition, where PSB typically occurs, this condition 


can be expressed as /ill > Ad ( 1,1) > / i r ( 0,2) > Hrl 
for transport in the forward (1,1) -A (0, 2) direction and 
Tll < Ti{ 1 , 1 ) < Ur{ 0,2) < /irl for transport in the re¬ 
verse (0,2) -A (1,1) direction. These conditions describe 
regions of transport known as bias triangles, with sizes 
proportional to the applied bias voltage, V^ as .[9] 

For a more complete picture of transport through a 
DQD near the (1,1) «A (0,2) charge degeneracy we 
now also consider the spin state of the electrons. For 
transport in the forward (left to right) direction cur¬ 
rent flows via tunneling events that cycle through the 
(0,1) -A (1,1) -A (0,2) -A (0,1) charge configurations. 
In the reverse direction it occurs via (0,1) -A (0,2) -A 
(1,1) -A (0,1). In both cases the first tunneling event 
brings an electron from a lead onto the DQD form¬ 
ing a two-electron system. These electrons can form 
one of four possible spin states, the s = 0 spin singlet 
| S) = (|tl) — lit))/A or one of the three s = 1 
triplet states |T 0 ) = (|TI) + |IT» /V2, |T+) = |tt), 

|T_) = |tt). In zero magnetic field, the | S) state 
is typically the ground state and the three degenerate 
triplet states are higher in energy by the exchange en¬ 
ergy J[10]. Since J scales with the wavefunction overlap 
between the two electrons, the charge configuration of the 
electrons will have a large influence on the magnitude of 
J. When in the (0,2) charge configuration, the overlap is 
large so J(0, 2) is as well, whereas in the (1,1) configura¬ 
tion the overlap is minimal as is J(l, 1) [11]. This differ¬ 
ence in exchange between the two charge configurations 
is the cornerstone of PSB. It leads to a situation where 
although the forward (1,1) -A (0, 2) interdot transition is 
allowed for the | S) ground states, current does not flow. 
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This occurs because, with J(l, 1) < J(0, 2), it is possible 
that even though conduction through the singlet states 
is allowed, i.e. [i LL > /a\ S \ 1,1) > //[^(0,2) > fi RL , 
the interdot transition for the triplet states is forbidden, 
//j Ta ^(l,l) </i[ Ta ^(0,2). Thus, once loaded, the |(1,1)7k) 
can tunnel neither forward nor back until it has relaxed 
into a singlet state. Since the spin relaxation from a \T X ) 
state to a | S') state is typically long relative to the tun¬ 
neling times, the measured current in the PSB region is 
near or below the noise floor of the measurement system. 
This situation truncates the size of the bias triangles in 
the forward bias direction, as shown in figure lc, while 
leaving the triangles measured in the reverse bias direc¬ 
tion at “full size”, as in figure lb. This size difference is 
the most basic signature of PSB and has been measured 
by several groups in various material systems[5, 12-21]. 


For a more intuitive presentation of the voltage space 
measurements depicted in figure 1(b) and (c) one can 
convert them to chemical potential space, or //-space [22], 
where triangles from both bias polarities can be shown on 
a single plot. Plotted in figure 1(d) is the //-space equiv¬ 
alent to what is shown in both panels (b) and (c), where 
the axes are defined as A/q(l, 1) = /i\ S \ 1,1) — Ull and 
A// r (0,2) = fi\?\ 0,2) — Ull- Also shown in this plot 
are outlines of the “state triangles”, which correspond to 
the conditions //ll > fif tate (l, 1) > fi s r tate { 0,2) > / lrl , 
in the forward biased direction and /i R l > fi^ tate ( 0, 2) > 
//State ( 1 , 1 ) > with /Ill = 0, in the reverse direc¬ 

tion, where state refers to the spin state, | S) or |T X ), of 
the two electron system. In figure 1(c) and (d), the re¬ 
gion labeled with the blue diamond corresponds to the 
region where PSB occurs, as shown in the corresponding 
chemical potential diagram. 


When a magnetic field 5, is applied, several effects 
can occur depending on the direction of the applied field. 
If applied perpendicular to the two-dimensional electron 
gas (2DEG) there can be a significant change in the spa¬ 
tial wavefunctions of the electrons. This will in turn af¬ 
fect the magnitudes of 7(1,1) and 7(0,2) leading to a 
change in size of the PSB region[12]. At higher fields, 
the | T x ) states can fall below the \S) states leading to 
PSB in the (1,1) —>• (0,2) direction caused by singlet 
states rather than the triplet states [23]. When B is ap¬ 
plied in the plane of the 2DEG the spatial effect on the 
wave functions is negligible and observed effects are due 
to the Zeeman splitting of the triplet levels. This in¬ 
plane magnetic field dependence will be the focus of the 
remainder of this manuscript. 


When considering the in-plane 5, we make the follow¬ 
ing assumptions for 5-dependence of the chemical poten- 



FIG. 1. Traditional forward Pauli-spin blockade, B = 

0 T Bias triangles at the (1,1) (0,2) charge transition 

with no applied magnetic field, Vbias = ±1 mV and J(0, 2) = 
0.3 meV. (a) Schematic of a typical double quantum dot. 
A bias applied between the left and right leads can lead to 
electrons tunneling through the device. The potentials of the 
dots are controlled with two capacitively coupled gates, (b) 
and (c) are voltage space pictures of current measurements 
in the (0,2) —y (1,1) and (1,1) —y (0,2) biasing directions 
using typical voltage and capacitance values. Filled regions 
of color correspond to non-zero current, (d) is the chemical 
potential space picture for the same transitions. Triangles 
below the A//j(l,l) s= 0 axis are for (1,1) —y (0,2) biasing 
and those above are for (0,2) —>• (1,1) biasing. Chemical 
potential diagrams for certain regions of interest in (b), (c) 
and (d) are also shown. The //-space origin is determined by 
the Fermi level of the left lead which we take to be at the 
ground potential for all measurements. The PSB region is 
labelled with the blue diamond. 
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When B is increased from zero, the first notable ef¬ 
fect is the change in size of the bias triangles. As re¬ 
ported in reference [13], the triangles in the forward 

(1,1) (0, 2) biasing direction grow in size proportional 

to the Zeeman energy Ez = ggsB, with g the electron g- 
factor, and gs the Bohr magneton. However, we predict 
a concomitant reduction in size of the triangles measured 
in the reverse (0, 2) (1,1) direction. Both of these ef¬ 

fects are illustrated in figure 2. 

Figure 2(a) shows the g -space picture of the chang¬ 
ing triangles (see supplemental material for details on 
using the g -space picture to determine the regions of al¬ 
lowed current). The origin is defined by the point where 
g\ S \ 1,1) = g[ S \ 0,2) = g LL = 0. Since g\ s \n : m) has 
no magnetic field dependence, our definition of the origin 
relates it to a single point in voltage space, at all mag¬ 
netic fields. This has the advantage that the direction 
the growth/contraction occurs along is clarified. The tri¬ 
angles for both biasing directions change size along the 
A/q(l, 1) axis. This is understood by considering the pro¬ 
cess which determines the edge of the triangles in each 
case. For the forward (1,1) —> (0,2) direction, the cru¬ 
cial process is the loading of an electron from the left 
lead onto the left dot. As B is increased, the | (1,1)T_) 
level is lowered by the Zeeman energy meaning that the 
loading of that state can occur at larger A/q(l,l) val¬ 
ues. This is shown in the chemical potential diagrams 
at the top of figure 2a for the point labelled with the 
blue triangle. In this region, since the only accessible 

(1,1) state is |T_), conduction is spin-polarized. In the 
reverse (0,2) —>> (1,1) direction, the crucial process is 
the Coulomb blockade of current due to loading of the 
|(1,1)T_) level below gLL- Again, as B is increased this 
level drops below gLL at a higher A/q(l, 1). This is de¬ 
picted in the chemical potential diagrams at the bottom 
of figure 2a for the point labeled with the red square. In 
both biasing directions, when viewed in this //-space pic¬ 
ture, the direction the triangles are changing is clear, and 
we can then easily generate the correct gate voltage space 
picture. This is done in figure 2(b) through (d), using 
typical capacitance and voltage parameters, which show 
bias triangles for the forward (1,1) —>> (0,2) direction 
growing with magnetic field in roughly the — Vi g direction 
(as compared to the detuning direction i.e. perpendicu¬ 
lar to the base of the triangles). This prediction is clearly 
observed in reference [13] and validates the assumptions 
made in equations 1 through 3. 

As B is increased further there will be a crossover from 
PSB in the forward (1,1) —)► (0, 2) direction, as discussed 
above and in previous works[1, 5, 12-21, 23, 24], to PSB 
in the reverse (0,2) —>• (1,1) direction. The onset of 
this reverse PSB is accompanied by the development of a 
tail-like feature in the regions of allowed current [25] and 
occurs at a magnetic field B$t where Ez = J(0, 2). 

Both of the high field features can be understood by 
considering figure 3, a //-space picture of a transport mea- 



FIG. 2. Magnetic field dependence for B < Bst: (a) g- 

space picture of bias triangles in both biasing directions with 
Vbias = ±1 mV, J(0, 2) — 0.3 meV at three B values corre¬ 
sponding to Ez — 0.0, 0.125 and 0.25 meV respectively. The 
triangles outlined with consecutively thicker lines correspond 
to higher B. As B is increased, the triangles in the forward 

(1,1) —>• (0,2) direction (left set of triangles) grow along the 
gi( 1,1) direction while the reverse (0,2) —(1,1) biased tri¬ 
angles shrink in the same direction. The change in triangle 
size is due to the Zeeman splitting of the | (1,1)T_) level as 
shown in the chemical potential diagrams corresponding to 
the points labeled with the blue triangle and the red square 
respectively, (b) through (d) are gate voltage space bias tri¬ 
angles in the forward (1,1) —» (0,2) biasing direction at the 
same magnetic fields as (a). The vertex labeled with the star 
occurs at the same position at all magnetic fields in both g- 
space and gate voltage space. The growth direction we predict 
is clear and agrees with the data in Figure 4 of reference [13]. 


surement in the reverse (0,2) —>> (1,1) direction with 
Ez > J{ 0,2). Included are the state triangles for the 
various singlet and triplet chemical potential levels. The 
regions labeled with the blue diamond and square corre¬ 
spond to the region of PSB. The physics of this PSB is 
similar to the traditional low field (1,1) —>> (0,2) PSB, 
but with the | (0, 2)5) state filling the role of the trapped 
state and the |T_) levels contributing to current. This 
is shown in the corresponding chemical potential di¬ 
agrams. Throughout the PSB region, although cur¬ 
rent is allowed through the ground state |T_) levels, 
the |(0,2)5) state can load from the right lead and, 
since the |(0, 2)5)—»|(1,1)T_) tunneling event is not al¬ 
lowed, it will block current until it relaxes to a |T_) 
state. Much like how the leakage current in the low 
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FIG. 3. B > Bst //-space picture of transport in the re¬ 
verse (0,2) —> (1,1) direction with an applied magnetic field 
gUBB — 0.5 meV, J(0,2) — 0.3 meV and Vbias — —1 mV. 
The region of allowed current (red) is no longer triangular. 
The tail-like feature, indicated with the star, consists of cur¬ 
rent through only the |T_) levels and therefore is a spin po¬ 
larized current upon exiting the DQD. The regions identified 
with the diamond and square symbols define the area of PSB. 
Unlike the typical PSB, this high-field PSB is caused by the 
loading of the 1(0,2)$) level rather than the |(1,1)T) level. 
This means that any leakage current measured in this region 
is due to relaxing from the |(0, 2)$) and |(0, 2)T_) states. In 
the region identified by the black circle current is allowed 
through both the | S) and |T_) states. Below the //-space plot 
are chemical potential diagrams for the indicated positions 
with key transitions indicated with arrows. 


field (1,1) —> (0,2) PSB region was used to probe the 
|(1,1)T)-H(1,1)S) relaxation[5, 13, 14, 16, 18, 21, 24], 
leakage currents in this reverse PSB region can be used 
to probe the |(0, 2)S)—»|(0, 2)T_) relaxation. 

Figure 3 has another distinct feature, the tail region 
labeled with the star. This feature also develops when 
the |T_) state becomes the ground state in the (0,2) con¬ 
figuration. The feature can be understood by considering 
the border between the tail and the PSB region. To the 
left of this border no current is measured because the 
|(0, 2 )S) can load and prevent current flow as previously 
discussed. To the right of the border the | (0, 2 )S) level is 
above /irl and is therefore unable to load (compare the 
chemical potential diagrams labelled with the diamond 
and the star). Since the |(0, 2)T_) state is the only state 
that can load, the current measured in the tail region will 
be spin polarized upon exiting the DQD, which suggests 
a potential use as a polarized spin current source. [26] 

By thoroughly considering the //-space picture typi¬ 
cally used to explain PSB, we have identified several in¬ 


teresting and unexplored features involved in the (1,1) 

(0, 2) transition of a double quantum dot at both low 
and high magnetic fields. The boundary between the 
two field regimes occurs at Bst where the Zeeman en¬ 
ergy of |(0, 2)T_) equals J(0, 2). With g = 2 and ex¬ 
change energies reported[13, 17, 20] in the range of 0.24 
to 1.4 meV, Bst in Si should be in the range of 2 to 
12 Tesla, accesible with standard cryomagnetics. This 
suggests that Si offers an opportunity to test our predic¬ 
tions. Simple validation can be achieved by performing 
measurements of the type reported in [13] for both bi¬ 
asing directions and comparing with our predictions of 
changing triangle size. Following that, the high-field tail 
feature we predict should be a clear signature that the 
reverse PSB regime has been reached. This naturally 
leads to experiments involving the leakage current in the 
PSB regime. Similar to low-field PSB leakage current 
investigations^, 5, 13, 14, 16, 18, 21, 24], these mea¬ 
surements would probe singlet-triplet relaxation mech¬ 
anisms but in the (0,2) charge configuration rather than 
the (1,1). This removes any spin-orbit interaction and 
nuclear field gradients between the two electrons, both of 
which were previously shown as the dominant coupling 
mechanisms between the | S) and | T x ) states in certain 
material systems[5, 14, 16, 18, 21, 24, 27]. Thus, mea¬ 
surements of this type will deepen our understanding of 
higher order spin relaxation mechanisms. 

We are grateful for the many useful discussions with 
Josh Pomeroy, Garnett Bryant, Mark Stiles and Michael 
Gullans. 
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( 1 , 1 )->( 0 . 2 ) 



FIG. SI. Low field state triangles: State triangles in the 
(1,1) —>> (0,2) direction in //-space with J(0, 2) = 0.3 meV, 
Ez — 0.15 meV, and Vbias — — 1 mV. The singlet trian¬ 

gle always lies at the origin by definition, and all other state 
triangles are shifted along the two axis according to the ap¬ 
propriate energy splittings. 

ALLOWED CURRENT WHEN J(0, 2) > E z 

In figure 2 of the main text we show the region of 
allowed current in a //-space picture for low magnetic 
fields. Here, we elaborate on the details of where and 
why current is allowed. The process of making chemical 
potential space plots is rather straightforward. For each 
state, one right-angle triangle is formed with horizontal 
and vertical sides of length eV^ias • The singlet triangle 
has the appropriate vertex positioned at the origin while 
the other state triangles are shifted according to their ap¬ 
propriate exchange and Zeeman energies. This is shown 


in figure SI with Vbias — —1 mV, J(0, 2) = 0.3 meV, and 
Ez = 0.15 meV. 

To determine the regions of allowed current one can 
draw the state triangles for the magnetic field of interest 
and then examine the different regions they define. By 
counting the conduction paths formed by transitions that 
are both energetically favorable and allowed through spin 
conservation one can determine how many channels are 
allowed in a given region. Secondly, one must determine 
if there are any partial channels, that is states that can 
be loaded onto the double dot but not contribute to fur¬ 
ther tunneling, such as the |(1,1)T X ) levels in traditional 
PSB. If one of these states exists current is not expected. 
Figure S2 shows this process for the //-space plot in fig¬ 
ure SI. The filled regions are colored according to the 
number of channels conduction can occur through. The 
red region is the PSB region, where, although conduc¬ 
tion is allowed through the singlet channel, no current 
is expected to flow due to the loading of the | (1,1)T_) 
state. In the green region all four states contribute to 
conduction. The cyan region corresponds to where the 
|(1,1)T + ) level cannot be loaded but all other states con¬ 
tribute to current. Finally, the yellow region corresponds 
to where there are no blocking partial channels and only 
the |T_) levels contribute to conduction. 


INTERNAL STRUCTURE 

Looking at figure SI, one might expect some current 
dependence relating to the number of conduction chan¬ 
nels allowed. Any dependence would likely depend heav¬ 
ily on any state dependent tunnel rates and very device 
specific. 
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FIG. S2. Low field state triangles: State triangles in the 
(1,1) —> (0,2) direction in /x-space with J(0, 2) = 0.3 meV, 
Ez = 0.15 meV, and Vbias — — 1 mV. Colors in the filled re¬ 
gion correspond the number of allowed conduction channels. 
The red region is the Pauli-spin blockade region with one al¬ 
lowed channel but one blocking partial channel. In the green 
region conduction is allowed through all four states, | S), |To), 

| T+) and |T_). In the cyan region the |T+) level cannot load 
so it is not involved in conduction whereas in the yellow region 
only the |T_) level can load. Thus, current is polarized. Also 
shown are chemical potential diagrams corresponding to var¬ 
ious points in the plot, which indicate whether or not current 
is allowed. 
















































